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Abstract 

During heavy ion collisions, high temperatures and strong magnetic fields are generated. We employ the gauge- 
gravity duality to study the Nf = 2 QCD phase diagram under these extreme conditions in the quenched appro- 
ximation, in particular we use the non-antipodal Sakai-Sugimoto model (SSM). We take the different coupling 
of up and down flavours to the magnetic field into account geometrically, resulting in a split of the chiral phase 
transition according to flavour. We discuss the influence of the magnetic field on the chiral temperatures — 
in physical GeV units — in terms of the choice of the confinement scale in the model, extending hereby our 
elsewhere presented discussion of fixing the non-antipodal SSM parameters to the deconfinement phase. The 
flavour-dependent (T,L,eB) phase diagram, with variable asymptotic brane-antibrane separation L, is also pre- 
sented, as a direct generalization of the known (T,L) phase diagram of the non-antipodal SSM at zero magnetic 
field. In particular, for sufficiently small L we are probing a NJL-like boundary field theory in which case we 
do find results very reminiscent of the predictions in NJL models. 



1 Motivation 



It is already known for some time that gigantic magnetic fields occurred during the cosmological electroweak 
phase transition [lj, but recently, it has also become clear that magnetic fields up to 10 15 Tesla can occur shortly 
during relativistic heavy ion collisions, see e.g. ElO. To put such a field strength in the correct perspective, a 
magnetar, or highly magnetized neutron star, reaches "merely" 10 9 Tesla. The interest in QCD studies under these 
extreme conditions has therefore increased considerably. As strong coupling effects are relevant in the setting 
of interest, we must rely on nonperturbative tools to study the relatively new realm of QCD physics in a strong 
magnetic background. This has initiated a vast amount of original research, let us refer to HI and references therein 
for a recent review. In particular, confinement and chiral symmetry breaking are two typical nonperturbative QCD 
effects that can be affected. It is common knowledge that QCD deconfines at a certain temperature T c , while 
chiral symmetry is restored at T%. We shall be concerned with QCD in the chiral limit here, i.e. we ignore the 
bare quark masses, in which case a clear-cut chiral transition exists. In real life QCD, with massive dynamical 
flavours, only approximate order parameters can be defined for both chiral restoration (~ chiral condensates) and 
deconfinement (~ Polyakov loop), leading to cross-over behaviour rather than sharp phase transitions. Once a hot 
debate whether T c and T % coincided or not, see JUO for 2 views on the Nf = 2 + 1 case, it is by now accepted they 
are close in the cross-over region. These results were obtained using lattice simulations, a powerful ally to access 
the nonperturbative QCD sector. 
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More recently, a vivid discussion evolved around the possibility that T c and T x separate under the influence 
of a constant magnetic background field, B = Be z , and this for Nf = 2 QCD, see e.g. [H |9). Although the 
Nf = 2 lattice results of [9] indicated a weak rise in the transition temperatures T c and T x , both remained com- 
patible with each other (a split of ~ 2%), while the various analytical model based results were inconclusive on 
the matter, as different results were obtained per QCD model (71 [H. Somewhat later, a more thorough lattice 
study appeared using Nf = 2+1 flavours with physical masses, leading to a much more complicated behaviour 
in the chiral/deconfinement (pseudo-)order parameters and ensuing critical temperatures JTOl . It was motivated 
that the reported behaviour — where contrasting with the results of 19] — should be traced back to the lighter dy- 
namical flavours and partially also to the present strange flavour, as the up (u) and down (d) quark of [9] were 
considerably heavier. Soon after, the first analytical papers appeared trying to explain the state-of-the art lattice 
data using backreacting pion dynamics ifTTTl . The naive reason for expecting a split between T c and T % was the ex- 
pected enhancement of chiral symmetry breaking due to a magnetic field — the so-called chiral magnetic catalysis 
|[T2l . The results of [10] showed a more subtle picture: the magnetic catalysis was confirmed for temperatures 
(sufficiently) below T c , but for larger temperatures the (averaged over up and down) chiral condensate displayed a 
non-monotonous shape, a feature translated into a similar behaviour in the transition temperature. This observation 
of an "inverse magnetic catalysis" seems to depend crucially on taking into account quark backreaction effects, 
see also [ 13], so we do not expect it to appear in the unquenched Sakai-Sugimoto setting we will use. 

Since a magnetic field couples to the up and down flavours with another strength, as they carry different electric 
charges, it seems natural that the up and down chiral restoration temperatures can be different, as well as the mag- 
netic catalysis itself. We recall that the classical chiral structure of QCD with and without magnetic field is diffe- 
rent, since coupling a magnetic field to the quarks reduces U (2) L x U (2) R to [E/(1)l x U(1)r] u x [E/(1)l x U(l)R] d , 
so that the eventually broken chiral invariances U(\) U A and U(\) d A can experience a different restoration tempera- 
ture. Lattice simulations indeed confirm a larger value for the (uu) than for the (dd) chiral condensate at T = 
|[T4l . as does the Nf = 2 Nambu-Jona-Lasinio (NJL) model JHIl . It would appear natural that should conse- 
quently be larger than The splitting of degenerate order parameters, like (uu) and (dd) at eB = 0, when an 
external field is switched on, is not that unfamiliar. In certain exotic superconductors, e.g. S^RuCU, a similar 
phenomenon occurs |[T6l . 

Here, we will use the gauge-gravity duality, a powerful tool to analytically study certain aspects of strongly 
coupled gauge theories JTTIl . to shed further light on the possible (T,eB) QCD phase diagram. In particular, we 
rely on the non-antipodal Sakai-Sugimoto model (SSM) llT8ll with 2 quenched flavours, represented by two D8-D8 
probe brane-antibrane pairs, in the chiral limit and with 3 colours. Going beyond the quenched approximation 
and/or including massive dynamical flavours leads to utter complications H9H20). We take into due account the 
different coupling of up and down flavour to the magnetic field, leading to a different probe brane geometry per 
flavour. We compare our findings with the lattice results of [22]. Although those results are referring to two-colour 
QCD, it contains an extrapolation to the chiral limit, which is the closest the available lattice results come to the 
Nf = 2 SSlvQ In ll22l . no manifest split between T c and T x was reported, whilst the chiral condensate increases 
monotonically with the applied magnetic field for all temperatures in the confinement phase. We do present similar 
results here using prefixed values for the string theory parameters of the SSM. With prefixed, we mean that a few 
physical QCD input values at zero magnetic field are chosen to match the corresponding SSM predictions. These 
results are a generalization to the 2 flavour case of the single flavour analysis of [23]. To extend the scope of 
our analysis, we will also allow that the asymptotic D8-D8 separation L — or equivalently, the confinement scale 
M — can vary and as such we construct the magnetic generalization of the (T,L) phase diagram of Figure 7 in the 
original work ll24l concerning the SSM phase diagram. All results are presented in GeV units to make comparison 
with other QCD approaches more direct. 

1 We will use N c = 3 colours to get explicit numbers, where in principle the limit N c — » °o is always understood at the holographic level. 
Setting N c = 2 or N c = 3 will not induce any qualitative change in the SSM results. 
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2 Setup 



2.1 The Sakai-Sugimoto model 

At zero temperature, the SSM llT8l involves a system of Nf pairs of D8-D8 flavour probe branes placed in the 
D4-brane background 

ds 1 = Q 3/ \ n ^dx v + f(u)dx 2 )+(*y 2 (J^+u 2 dCli\, /(„) = l-4 > (i) 

where R 3 = Kg s N c £ 3 , with g s (£ s ) the string constant (length). There is a natural cut-off at u = uk, which ensures 
confinement in the dual field theory living at the boundary u — » <*>. A smooth cut-off is realized if x has a periodicity 
Sx = ^R^^Uf}^ 2 = 27iM _1 with M the compactification scale. A U(Nf) L x U(N/)r gauge theory resides on the 
stack of coinciding D8-D8 flavour pairs, which corresponds to the global chiral symmetry in the dual QCD-like 
theory. The cigar-shape of the (w,x)-subspace of the D4-brane background forces the embedding of the flavour 
branes to be U-shaped, which signals the breaking of chiral symmetry U (A^)l x U (N/)r —> U (Nf)v as the merging 
of the D8-branes and D8-branes at the value u = uo > uk. The value uo is directly related to the asymptotic 
separation L (at u — >> ©o) between D8- and D8-branes, indicated in Figure [l] For values of uo greater than or equal 
to uk the embedding is respectively non-antipodal or antipodal. 



U/K 




Figure 1: The Sakai-Sugimoto model. 



At finite temperature there are two regular Euclidean backgrounds with the same asymptotic geometry that 
compete with each other in the partition function, the Wick-rotated version of the D4-brane background ([T]) and a 
black D4-brane background l24l . The Wick-rotated D4-brane background has a cigar-shaped (w,x)-subspace and 
a cilinder- shaped (w^)-subspace, with the periodicity of the /-circle arbitrary and equal to the inverse temperature, 

(3 = and the periodicity of the x-circle fixed to St = ^/? 3 / 2 w^ 2 , while the black D4-brane background 

ds 2 = (-) 3/2 ' (f(u)dt 2 + 8^(1x^x1 +dx 2 )+f-V (if^ + u 2 dQl), /(V) = l-^|, (2) 
\RJ \uj \f{u) J u 5 

has a cigar-shaped (w,/)-subspace and a cilinder- shaped (w,T)-subspace, now with the periodicity of the T-circle 

arbitrary, but the periodicity of the /-circle fixed to 8/ = (3 = T~ l = ^R 3 / 2 u^ 1 ^ 2 . These two backgrounds are 
identical, modulo a redefinition of coordinates / and x, when St equals (3, which happens at the deconfinement 
transition temperature T c = ^R~ 3 ^ 2 u)/ 2 . In the deconfining phase 2 ] the embedding of the flavour branes is no 

2 Let us remark here that in H25I some problems concerning the identification of the deconfined phase with the black D4-brane back- 
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longer forced to be U-shaped, as the (w,x)-space is no longer cigar-shaped. At a certain value of the temperature, 
T X >T C , it will become energetically favourable for the flavour branes to fall straight down instead of merging in 
a U-shape, indicating chiral symmetry restoration (see Figure [2]). The essential features of QCD, chiral symmetry 
breaking and confinement, plus the ensuing chiral restoration and deconfinement at sufficiently large temperatures 
are thus nicely resembled by the SSM, amongst other QCD phenomenology JT8ll . To make further explicit contact 
with QCD, we determine the string related parameters of the SSM in physical (GeV) units, so that an explicit 
comparison with other approaches comes within range. 




Figure 2: (a) Deconfinement transition at T c and (b) chiral symmetry restoration at T x (> T c ). 



2.2 Fixing holographic parameters at eB = 

In iTTSl the independent parameters M and K of the antipodal (uo = uk) SSM were fixed to M ^ 0.949 GeV and 
K = X/(72n 3 ) w 0.00745, with X = gj M N c the 't Hooft coupling, by matching to the QCD input values 

f n = 0.093 GeV and m p = 0.776 GeV (3) 

for the pion decay constant f % and the p meson mass m p . Due to the relations uk = M _1 , R 3 = \M~ 3 and 
£7^7 3 = ^jM 3 1 18] all other parameters of the model are then fixed as well. Note that throughout this paper we 

ground were discussed, and instead a different background was proposed, namely a localized D3-soliton geometry. To make calculations 
of the critical temperatures feasible, it is however necessary to consider a high-temperature approximation of that background. There are 
moreover some subtleties concerning the inclusion of flavour branes, but the end result for the (T,L) phase diagram, see Fig. 10 of ['25], 
gives a qualitatively similar result as the 'old' SSM. Based on this observation, we can expect that qualitative features of the ^-dependence 
of the chiral transition temperatures are not unlikely to be similar in both backgrounds. It would be interesting to check this explicitly, 
but for the reasons mentioned above we will consider the simpler black D4-brane background in this paper. We thank Takeshi Morita for 
discussion on this point. 
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setN r = 3. 



In the non-antipodal case (uo > uk) the matching conditions ([3]) do not fix all freedom, since there is one extra 
parameter present, uq. From the eigenvalue equation which determines m p holographically, 



dz {^/^ ^^ = -\u^ul l2 iF?ml^ V p (0) = 0, y p (±oo) = 0, 



(4) 



with 



y (z) = J 3 * : 8 5-3T, u\ = ul + mz 2 , (5) 

V U z( U z - u k)-( u 0- U U k) 

we extract the values of uo that, for a given M, lead to m p = 0.776 GeV. The resulting function uo(M) is plotted 
in Figure [3] for a range of M — the maximum value of M corresponding to the limiting case uo —> uk — 1 /M — 
alongside with the function L(M) for the corresponding asymptotic separation between branes and anti-branes, 
determined from 




(6) 



u*f-u«f 

with f(uo) denoted as /o- Next, demanding that the SSM-prediction for the pion decay constant 




A(M,a ,K) = W-KM 7 / 2 - 2/ dz -, , , - ^tt2 , (7) 



equals 0.093 GeV, leads to the function k(M) of allowed values for k, as plotted in Figure 3(c) The string tension 
(2710c') _1 = 871 2 M 2 k(M) is then also known as a function of M. 

The remaining freedom of choosing the mass scale M can be fixed for example by matching the SSM-prediction 
for the constituent quark mass, 



m^(M,w ,K) 87r 2 M 2 K / da . . 



(8) 

(Mm) 3 

to a phenomenologically reasonable value, as can be read off from Figure |4j 
In lf26ll we opted to reproduce m q = 0.310 GeV, leading to the following set of fixed holographic parameters: 

M w 0.7209 GeV, L w 1.574 GeV" 1 and k = w 0.006778. (9) 

21671^ 

With these values, the effective QCD string tension is computed to be a « 0.19 GeV 2 ll26l , in excellent agreement 
with the standard lattice value, a « 0.18-0.19 GeV 2 , extracted from ETIl . Said otherwise, we could equally well 
have selected the string tension as our third QCD input value instead of the constituent quark mass. 

Here, we will however opt to leave M variable, or equivalently L via Figure |3(b)[ with the eye on drawing the 
(T,L,eB) phase diagram later, and, more importantly, with the idea that the choice of M or L should be left free, 
as it determines the choice of holographic theory: L very small (~ St = ^ large ~ M small) corresponding to an 
NJL-type boundary field theory ||3TJ[30l[24l versus L = Sx/2 maximal (~ M maximal) corresponding to a maximal 
probing of the gluon background, i.e. the original antipodal SSM. 
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(a) (b) (c) 

Figure 3: (a) Values of uo (GeV -1 ) and (b) corresponding values of L (GeV -1 ) compatible with m p = 0.776 GeV. 
(c) Values of K compatible with m p = 0.776 GeV and f % = 0.093 GeV 




(a) 



Figure 4: SSM-prediction for m q (GeV) as a function of the confinement scale M (GeV), compatible with m p = 
0.776 GeV and f n = 0.093 GeV. 

2.3 Applying a magnetic field 

The next step is to introduce a magnetic background. To this purpose, we take a closer look at the U (Nf) gauge 
field A m (x^, u) (m — 0, 1, 2, 3, u) living on the D8-brane, with the non- Abelian DBI-actioiQ B2 

S = -2T % J d 4 xdu J 8 4 STr^-det[^ + (27ca / )iF mn ], (10) 

where the factor 2 in front makes sure that we integrate (in u) over both halves of the U-shaped D8-branes, = 
g s (u/R) 3 / 4 , 84 is the volume form of the unit 4-sphere in the background, T% = \/ ((27i) 8 ^) the D8-brane tension, 
STr the symmetrized trace defined as STr(Fi - -F n ) = ^Tr(permutations of F\ • • -F n ), the induced metric on 
the D8-branes, a! = l 2 s the string tension, and F mn the usual field strength; we use anti-Hermitian generators. As 
explained in IH1I261, one can slightly gauge the global U (Nf)v symmetry in the boundary field theory, i.e. make 
g(u —¥ 00) e U (Nf)y x-dependent, so that gd^g~ l = A^u 00) corresponds to adding a background U (Nf)v field. 
Working in the A u = gauge, the flavour gauge field can then be expanded as iTTSl 

A^,u) = A A/ (x ;U )+rest, (11) 

where "rest" refers to pions and vector mesons in the boundary QCD theory, which are irrelevant for the current 
purposes. 

3 As the 't Hooft coupling X is large, we ignore the Chern-Simons part of the action in the analysis, being a factor l/X smaller than the 
DBI-action. 
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An electromagnetic background field A e ™ can be switched on through 

A" = -ieQ em A; m = -ie ( 2/ Q 3 _° /3 ) A" = -fe (| + f ) A e ™ , (12) 

3 q 3 

with Q em the charge matrix for the up- and down-quarks. The choice A 6 ™ = xiBd^, or A 2 = x\eB, A 2 = A 2 /3, 
ensures the desired constant magnetic background B = Be^. The corresponding field strength tensor reads 

T K _ aA-zfi? *) = -,( I- I). us) 



3 l 



e£ J \ F d 



3 Results 



3.1 eB-dependent embedding in confinement phase 

For completeness, we briefly summarize here the discussion of the ^-dependent embedding of the N/ = 2 flavour 
branes in the confinement phase, i.e. the determination of u' — du/dx for each flavour, presented in more detail in 

ea. 



The induced metric on the flavour branes is given by 
dsm = gZd^dx" (m,« = ()••• 8) 



= (-) V ^v + ^-j _+(-) ^J^ + (-j ^ (14) 

where the metric components in flavour space are assumed to be different, to allow for a different response of both 
flavour branes to the magnetic field: 

with the w-coordinate appearing in g^ 8 , resp. g^ 8 following the up-, resp. down-brane, thus u E [wo, M ,°°] or w E 
[wo,j, °°] • Inserting this metric-ansatz into the DBI-action ( [TO] ), along with the background, gives S ccw / = S M + £^ 
with 



S £ = -r 8 ^ 4 W2j[ M ^" 4 y^(|) 3 + ^V^, A i = l + (27ax / ) 2 F?^ , (16) 

where the index £ refers to the up or down flavour, u' = du/di, % = f d 4 x and V4 is the volume of the unit 
4-sphere. The DBI-action for the two D8-branes thus reduces to the sum of two Abelian actions, which is the 

geometric translation of the explicit breaking of chiral symmetry U (2) a ^t/(l)^x[/(l)^. 
Suppressing the ^-index for the moment, u' can be determined from the conserved "Hamiltonians" 
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by assuming U-shaped embeddings, u' = at u = uo (with A(uo) and f(uo) denoted as Ao and fo, for each flavour). 
In the confinement phase, the ^-dependent action on each flavour brane is then given by 



poo 

Juo,e irfAe - u^foAoj 

and the asymptotic separations L = 2fdu/u' become (suppressing I for notational convenience) 

3/2 



if du(-} f- 1 . 



L conf (u ,eB / _ , .... , , , l( B 

Juo \uj y u 8 fA-ulfoA 

Jo 



Uq/oAo 

)A 

(19) 



3v^ V ^0 ^y/fA-foA&W 

To model magnetic catalysis, we keep L con f (wo, eB) fixed, for each flavour, at its flavour-independent value L for 
eB = 0. From the viewpoint of the boundary field theory, the flavour branes are extensive objects that stretch out 
infinitely far, viz. from u = <*> to u = uo into the higher-dimensional bulk space. As such, from this asymptotic 
perspective, it would appear that it would cost an infinite amount of energy to move the branes at u = <*>. Keeping L 
fixed thus seems to be a sensible boundary condition, one also used in e.g. [29, 30], to probe the effects of the bulk 



dynamics in the presence of the external field. The value of L, ranging from y to negligible w.r.t. y , determines 
how much of the bulk dynamics is probed, ranging respectively from all to none. 

So, as we keep L con f fixed to its value L at zero magnetic field, from here on also explicitly writing the M- 
dependence of L con f through R and 2710c', we can extract the ^-dependence of uo and of the corresponding con- 
stituent quark masses m q ^ u and m q ^, which correspond to the mass ^ of a string stretched between one of the 
flavour branes and the cut-off at uk = M~ l ll24l : 

L c i onf (u 0l eB,M)=L =^ u ,t(eB,L,M) =^ m q/ (eB,L,M). 



Using the correspondence between L and M as plotted in Figure 3(b) such that m p = 0.776 GeV, we obtain 
m q {eB,M) or m q (eB,L). This turns out to be a rising function of eB for all choices of M or L, modeling chi- 
ral magnetic catalysis, as already noted in ll23l for the single-flavour case. The magnetic field thus indeed boosts 
the chiral symmetry breaking, reflected in a stronger bending of the flavour branes, and the breaking is stronger for 
the up- than for the down-flavour, see Figure [5(a)| This was expected, since the up-quark couples twice as strong to 



the magnetic field. For the particular choice of parameters ([9]), m q is plotted in Figure 5(b) Both constituent masses 



show a quadratic dependence on eB for small magnetic fields, which then becomes nearly linear, and eventually 
thrives to saturation at very large eB. This linear behaviour is in relative agreement with recent lattice estimates of 
the chiral condensate^] in H4l[T0l|22]|. A linear behaviour, inspired by chiral perturbation theory, was fitted to the 
SU(2) lattice data for the chiral condensate in [33], while SU(3) data was fitted with a (eB) 3 / 2 behaviour in 041 : 
both these studies used 1 quenched flavour. Let us however point out that the plots of the lattice chiral condensates 
in all cited works do seem to display a nontrivial curvature for small values of eB. 



The stability of the flavour branes' embedding as shown in Figure [5(a)] is discussed in 1261 . At sufficiently low 
values of eB, it is stable, but at higher values of eB, we should not necessarily trust our findings, as an instability 
could occur in the p meson sector at eB ~ m 2 ~ 30 (w 0.6 GeV 2 ), driving a condensation ll35l . Evidently, this 
would lead to a different underlying action, describing the theory in the condensed phase. If this happens, also 
other models should take it into account. 



4 The concept of the chiral condensate is not well-defined in the SSM, we therefore use the constituent quark mass as a measure for the 
chiral symmetry breaking. 
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u ad(eB) 

u o (eB=0) 




m q (GeV) 




(a) 



eB (GeV 2 ) 



(b) 



Figure 5: (a) Up and down flavour brane embedding in the presence of a magnetic field, (b) corresponding con- 
stituent quark masses m q ^ u (blue) and m q ^ (red) for parameter choice ([9]). 

3.2 ^-dependent embedding in deconfinement phase 



Next, we turn to the finite temperature case. The background at T < T c is identical to the zero temperature 
background up to the period (3 = T~ l of Euclidean time, so nothing changes as compared to the T = case. 
Things get more interesting once we enter the deconfinement region. We again have an induced metric on each 
flavour brane, 



RJ 



(fdt 1 + bijdx l dx J ) + - u l d£l\ + 



R\ 



3/2 



with periodicity of the /-circle given by 



u J 



& = r- 1 = ^*vV /2 . 

3 T 



K\3/2 

r) 



1 /M\-3 1 



f\R 



du 



(20) 



(21) 



from which one can determine the action in the deconfined phase, completely analogous to the derivation of the 
action in the confined phase. For temperatures T < 7^, the ^-brane's embedding remains U-shaped, with action 



\ 



i 



(22) 



My)y 3 A £ 



where c = -2T%%V^g s l R 3/2 , y = u/u ,e, yr/ = u T /m& My) = 1 - (yr^/y) 3 and f j = 1 -y\ r If T > T x j, the 
^-branes are falling straight down, u' = oo ? with action 



S T £ >T%/ = couJJ / dy yy/y 3 At. 
The chiral transition temperature T x j is the temperature for which ASt becomes zero [24 j, with 



AS(u ,eB,y T ) 



actionu_ s hape - action strai ght . 



(23) 



(24) 



The correspondence between uo and L in the deconfined phase is modified into (again suppressing the flavour index 
here) 



L dec (u ,eB,y T ) 



2fl 3 /2 



3 Juq 



1 J-l/2^1/2 



(25) 
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As before, we will hold the asymptotic separation fixed at its starting value L at eB = and T = 0. This allows us 
to determine the eB- and T-dependence of uo (from here on also explicitly writing the dependence on M, through 
R and 2na'): 

L dec (u^eB,y T ,M)=L => uo(eB,y T ,L,M). (26) 



From Figure 6(a) it can be seen that the one-to-one correspondence between uo and L con f is not preserved in the 
deconfinement phase, where each value of L dec corresponds to two possible values of uo, as long as it does not 
exceed its maximum possible value (i.e. as long as T < T % ). We numerically verified that the energetically favoured 
solution for uo is the largest one, consistent with the intuition that the lower-wo solution contains more energy as it 
probes a larger portion of the background. Keeping L fixed during the deconfinement transition causes a jump in 
uq, as well as in the constituent quark and meson masses ll29l . 





(b) 

Figure 6: (a) L con f (blue) and L dec (GeV" 1 ) for T = T c (red) and increasing values of T(> T c ) at eB = 0, (b) 
AS/u /2 as a function of y T for eB = (blue), 0.5 (purple) and 1.6 GeV 2 (yellow). Both figures for M = 0.7209 
GeV. 



With the expression found for uq, the expression for AS at fixed L is also known, 

AS(u (eB,y T ,L,M),eB,y T ,M) = AS(eB, y T ,L,M), 



(27) 



so the chiral temperature can be determined from the point where the U-shaped embedding breaks into separated 
branes, i.e. when AS = (see Figure [6(b)] ): 

AS(eB,y T ,L,M) =0 =^ y\{eB,L,M). 

The corresponding value of uq at the chiral transition is then given by 

= uo(eB,y^(eB,L,M),L,M) = ufa{eB,L,M). 

Plugging the obtained y\(eB,L,M) and Wq(^5,L,M) into the definition for the chiral temperature 



471 V # 3 471 



R 3 4ti 



y\{eB,L,MY 



lu$(eB,L,M) 



-^^ = T x (eB,L,M), 



we obtain T x (eB,L,M). 



From the parameter discussion at eB = we know the value of the fixed asymptotic separation given a value 
for M such that m p = 0.776 GeV (Figure 3(b)), hence we obtain T x (eB,M) or T x (eB,L), to be compared with the 
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deconfinement temperature T c —M/ (27l). The deconfinement temperature T c will not change as it is determined 
from the background D4-brane metric, which could only become ^-dependent when the backreaction of the D8- 
branes would be taken into account. Or, field theoretically: T c is ^-independent in a quenched setup, because the 
magnetic field can only couple to the neutral gluons indirectly via the quark interactions. For every choice of M (or 
L), T x (eB) rises with eB ("chiral magnetic catalysis"). But depending on the choice, there will or will not arise a 
split between T % and T c . Doing this for each flavour, we find T^(eB,L) and T^(eB,L), with consistently higher 
than TJ* for a given value of L, as expected. This leads to an intermediate phase where the chiral symmetry for up 
quarks is still broken while the chiral symmetry for down quarks is already restored: 

u(iy v x u(i) d v % u(\y v x (u(i) L x u(i) R ) d T i (U(l) L x U(1) R ) U x (U(1) L x U(l) R ) d , 

as sketched in Figure |7J 
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Figure 7: Embeddings in the deconfined phase with magnetic field. 

In Figure[8]the (T,M,eB) and (T,L,eB) phase diagrams of the two-flavour non-antipodal SSM are plotted. This 
generalizes the Nf = 1 SSM phase diagram in Figure 7 of ll24l to the Nf = 2 magnetic case. For set-ups with 
large values of M, namely M > 0.767 GeV corresponding to m q (eB = 0) < 0.274 GeV or L > 1.681 GeV -1 , 
there is no split between T x (eB,M) and T c = M/(2n), no matter how large the applied magnetic field is. This 
is a consequence of the saturation of the rising of T x with eB. In a SSM with M < 0.657 GeV, corresponding to 
m q (eB = 0) > 0.353 GeV orL < 1.473 GeV -1 , there is already a split between chiral and deconfinement transition 
before the magnetic field is turned on: T x (eB = 0,M) > T c , which becomes larger as eB increases. This regime is 
probably the least physically relevant, as the values for constituent quark masses are too large and the values for 
the deconfinement temperature smaller than 0.105 GeV, which is rather small compared to the chiral limit value 
we can extrapolate "by hand" from ll28l giving T c ^0.150 GeV. The third possible case is that the value of M 
is such, 0.657 GeV < M < 0.767 GeV (~ 0.274 GeV < m q (eB = 0) < 0.353 GeV or 1.473 GeV -1 < L <_L.681 

For 



10 



GeV -1 ), that T % (eB = 0,M) = T c but a split between T % and T c arises at some value of eB, plotted in Figure 
each of the above possible cases, an exemplary cross section of the (r,M, eB) phase diagram is shown in Figure[9j 
the middle one corresponding to the best matching parameters for reproducing a reasonable m q (eB = 0) « 0.310 
GeV, although the corresponding value for T c « 0. 1 15 GeV is still on the small side. A deconfinement temperature 
T c 0. 150 GeV would correspond to M = 0.942 GeV, very close to the antipodal value so in the regime where no 
split arises between T % and T c , leading however to an unphysicall}^] small m q [eB = 0) « 0.046 GeV. 

The papers (23l |30) also pointed out the splitting of the critical temperatures for the one flavour version of the 
SSM, but leaving the parameters of the SSM undetermined, making an explicit comparison with other approaches 
harsh. The explicit breaking of the global flavour symmetry by the different electromagnetic coupling of the up and 



5 This might be related to the shortcoming of the SSM (in the used form, not considering possible modifications as in 1 20 1) that the bare 
quark masses always remain zero. 
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Figure 8: (a) (GeV) (upper blue surface) and Tj* (GeV) (lower blue surface) as functions of eB (GeV 2 ) and 
M (GeV) compared to T C {M) (GeV) (purple), (b) same with M-dependence replaced by L-dependence compatible 
with m p 0.776 GeV 
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Figure 9: Cross sections of Figure[8]for (a) M = 0.65 GeV, (b) M = 0.7209 GeV and (c) M = 0.77 GeV, respectively 
corresponding to m q {eB = 0) = 0.357,0.310 and 0.272 GeV and T c = 0.103,0.115 and 0.123 GeV The appearance 
of a split between T x (GeV) (blue) and T c (GeV) (purple) depends on the choice of M, or equivalently L. 



down flavour is also now taken into account for the first time, leading to a split between the two chiral transitions 
themselves. 

The fact that a split between T % and T c can emerge only for sufficiently small values of the asymptotic brane 
separation L, i.e. sufficiently close to an NJL effective description of QCD, seems to be supported by NJL model 
calculations |8|, see also the discussions in El \15\ [36), that seemingly contrast with lattice data. Selecting the 
holographic parameters in a way that brings the SSM as close as possible to (the chiral limit of) QCD, rather 
leads to a picture of the form of Figure |9(b)| or |9(c)[ namely no split at all or a small split that only emerges at 



rather large values of eB. Our findings are in this perspective consistent with lattice data of ll22l . where a split was 
neither reported. However, we must also repeat here that our results are obtained in a quenched framework, hence 
important QCD effects at the level of transitions can be missing (e.g. pion loop effects). In particular, magnetic 
effects on the deconfinement temperature cannot be taken into account in the SSM without including backreaction 
of the probe branes on the background D4 metric. 



3.2.1 Remark on the antipodal SSM 



In the original antipodal Sakai-Sugimoto model, with uo = uk and the asymptotic separation L taking its maximum 
possible value, the embedding of the flavour branes is unaffected by the presence of the magnetic field. From 
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this we can conclude that the antipodal Sakai-Sugimoto model is unable to capture the magnetically induced 
explicit breaking of chiral symmetry, as well as the chiral magnetic catalysis. Chiral symmetry restoration and 
deconfinement coincide for all values of the magnetic field. 



eB 




Figure 10: Value of eB where (GeV) (blue) resp. (GeV) (purple) becomes larger than T c (GeV) for confine- 
ment scale values 0.657 GeV < M < 0.767 GeV. 



4 Summary 

In conclusion, we have investigated the phase diagram of the two flavour version of the non-antipodal Sakai- 
Sugimoto model in the presence of a temperature T and external magnetic field eB. In particular we payed attention 
to fixing the holographic parameters, presenting a discussion of how they can be fixed by matching to carefully 
chosen QCD input parameters, in order to be able to present the phase diagram and related results in physical GeV 
units. This makes comparison to other approaches more direct. We indeed could compare our results with lattice 
and NJL results, the main conclusion being that the SSM results are consistent with other quenched settings that 
are able to model chiral magnetic catalysis. 

The main results are presented in the (T,L,eB) phase diagram in Figure [8] and cross sections of that plot for 
different values of L in Figure [9] Here, L is the asymptotic separation between the flavour probe branes. Keeping 
L fixed serves as a boundary condition for the bulk dynamics, the effective boundary theory ranging from the 
NJL-type for small L to a chiral QCD-like theory where gluon dynamics are fully taken into account for maximal 
L. The value of L, i.e. the choice of the type of boundary model in a sense, determines if a split between the chiral 
and deconfinement temperature may or may not arise, as summarized in Figure [9] Due to the different coupling to 
the magnetic background of the u and d flavour brane, we also find a split between the separate chiral transition 
temperatures. 

It remains a challenge to construct a holographic dual of realistic QCD that could also describe the complicated 
finite temperature (above and below T c ) behaviour of the chiral condensate as found in the latest lattice results 
fTTOll . This would require taking backreaction of the flavour branes on the background metric into account. A 
holographic model that could be interesting to look at from this point of view, is the Kuperstein-Sonnenschein 
model ll38L in which chiral magnetic catalysis was observed L21J in the non-backreacted case, and which was 
extended to include backreaction (be it perturbatively and for vanishing magnetic field) in P9lK using a smearing 
technique. A downside of this model is however that it does not incorporate confinement due to the choice of the 
background in which the flavour branes are placed. 

We end by noticing that the magnetic phase diagram could become even more intricate, as demonstrated in |[35l : 
the QCD vacuum could become superconducting at sufficiently strong magnetic field. This matter is currently 
under investigation also in the holographic framework, see also ESI [37). We hope the presented results will 
stimulate further research in the area of QCD in a strong magnetic background, where many exciting new physics 
could be awaiting discovery. The gauge-gravity correspondence can offer viable input, next to QCD model and 
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lattice based approaches. 
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